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BIO1819P0090US 

^ ^^^^^ MICROSCOPE SYSTEM FOR READ-TIME SIMULTANEOUS TPMPOu a t 
MEASUREMENTS OF METACHRONAL WAVE PERIO^S^^S^y B^?^ 

FREQUENCY 

FIELD OF THE INVENTION 

This invention relates generally to diagnostic systems and, in particular, to a method and 
apparatus for measuring ciUary beat frequency and metachional wave frequency in ciliated 
native tissues, ciliated organ cultures, ciliated monolayer cultures, and air-liquid interface ciliated 
cultures derived from vertebrates and invertebrates. 

GOVERNMENT SUPPORT 

The present invention was made with U.S. Government siqjport fiom the National 
Institutes of Health, undo: Gnmt No. HL 59126. The U.S. Government has certain rights in this 
invention. 

3ACKCa?.OUND OF THE INVENTION 

The effectiveness whereby ciliary activity of the ciliated cells in the airway epithelium 
adequately transports mucus through the airways is governed by both the fiequency of the cilia 
beat (ciliary beat frequency. CBF) and tfie efficiency of the antiplectic ciliary coordination as 
indicated by the metachronal wave (metachronal wave frequency. MWF). Conceptually, the 
mucociliary wave is composed of a low frequency signal, the metachronal wave, that is carried 
by a high frequency signal, the ciKary beat frequency [Sanderson, MJ, and MA Sleigh. Ciliary 
activity of cultured rabbit trachea cUiary beat pattern and metachiony. J Cell Sci. 47: 331-347, 
1981.]. Ihe MWF is typically determined from the elapsed time of the horizontal propagation of 
the ciliary wave between defined vectorial positions [Wong, L.B., Miller, I.R and Yeates, D.B. 
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The nature of mammaUan ciliary metachronal wave. J. Appl. Physiol.. 75(l):458-467. 1993], 
whereas measurement of CBF is based on the elapsed time it takes for the ciUa to return back to 
their starting vertical position [Wong, L.B.. Miller, I.F1 and Yeates. D.B. Regulation of ciliary 
beat frequency by autonomic mechanisms: in vitro. J. Appl. Physiol., 65(4): 1895-1901, 1988.]. 
Based on these principles, various techniques have been developed to either analyze pre- 
recorded images of ciliary motion off-line using statistical methods [Hennessy, S. J.. L.B. Wong, 
D.B. Yeates, and I.F. Miller. Automated measurement of ciliary beat frequency. J. Appl. Physiol. 
60(6): 2109-2113, 1986; Ben-Shimol, Y., L Dinstein. A. Mdsels and Z. Priel. Ciliary motion 
features from digitized video photography J. Compu. Assisted Microsc. 3:103-116, 1991.] or by 
usini dual focal spot laser light scattering techniques to perform real time analysis of the phase 
shift of the scattered signals [Wong, L.B., MiUer, I.F. and Yeates, D.B. The nature of 
mammalian cihary metachronal wave. J. Appl. Physiol., 75(l):458-467, 1993.]. These 
technologies are labor and time intensive, and do not necessarily analyze the CBF and MWF 
from the same group of dlia. If accurate analysis of CBF, MWF, and their coupling is to be 
achieved, it is necessary to develop a method to perform such an analysis on data from a single 
group of ciUa. Such a method would facihtate the determination of how the ciUa generate a 
horizontal prcipagating "staiiding" metachronal wave from a nonstationary, non-propagating 
asymmetrical vertical oscillatory motion, the CBF. The invention described herein allows such 
an analysis to be performed. 

With the assumption that the temporal signal derived from the ciliary activity at a single 
location is composed of both CBF. and MWF, one can mathematically derive both CBF and 
MWF if the spatial coherency of the focal spot has been assured. However, to resolve likely 
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nonstationary multiple frequencies from a single source in real-time without prior information of 
the nature of the signals, efficient algorithms need to be developed. 

Accordingly, there is a need for a new method and system which measures ciliary activity 
and that enables continuous real-time analysis of both ciliary beat frequency and metachronal 
wave frequency m excised ciUated tracheal tissues as well as in primary and secondary ciliated 
cultures. Such a method and system, described herein, utiUzes the concept of time-scale wavelet 
analysis and Hilbert transformation for backscattered light derived from a confocal (conjugate) 
spot on the moving ciliated surface. This light has inherent high and low frequency components 
corresponding to CBF and MWF. 

SUMMARY OF THE INVENTION 

The present invention relates to a method and system that enables continuous real-time 
analysis of both ciliary beat frequency and metachronal wave frequency from a single spot in 
excised native ciliated epithelial tissues as well as in primary and subsequent epithelial cultures. 
Such method and system utilizes the concept of time-scale wavelet analysis and Hilbert 
Transformation for backscattered ligjit derived from a confocal (conjugate) spot on the moving 
cilia. Tins light contains inherent high and low frequency component? corresponding to CBF 
andMWF. 

An object of the present invention is to enable remote detection and analysis of ciliary 
beat frequency and metachronal wave frequency from a single focal spot on either native 
epithelia, air-liquid interface cultures, organ cultures, or single cell culture. This is accomplished 
by designing a microscope-based system that enables the analysis of detected signals from the 
single focal spot in a way that preserves spatial coherency of the dynamic moving cilia. The 
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focal' spot of the microscope objective is the area of the primary Airy Disc. A pinhole 
conesponding to the area of the primary Airy Disc strategically placed in the detection and 
excitation optical paths serves as a conjugate (confocal) spatial filter to ensure that the detected 
signal from the focal spot is spatially coherent. 

The preferable configuration of the invention has an excitation and detection system 
comprised of 1) a monochrbmatic visible wavelength laser light source of beam diameter greater 
than 4 mm to enhance focal depth of the objective; and 2) a high speed photon detector to collect 
the backscattered light containing the signal. This is accomplished by splitting an expanded 
incident laser beam using a beam splitter with a minimum Lambda/4 unifonn flatness. The 
reflected laser beam is directed to the microscope objective. To provide the reference signals for 
hetCTodyne modulation, the transmitted laser beam i^ incident on a random depolarizer. 

Real-time analysis of the simultaneous nonstationary CBF and MWF signals is achieved 
using a computer controlled system with the appropriate software to analyze the photon detector 
outputs. It has been determined fiofti the piibKshed literature using laser light scattering 
techniques to determine, CBF, that the optimal sampling time is 2 to 6 milliseconds for 512 
channels [Chandra, T., Yeates, D.B., Miller, I J. and Wong. L.B. Stationary and nonstationary 
correlation-firequ^cy analysis of het«rodyne mode laser light scattering: Magnitude and 
periodicity of canine tracheal ciliary beat fi»quency in wvo. Biophys. J., 66:S78-890, 1994.], 
i.e., 1.5 seconds of data collection time is needed to guarantee the robust analysis of CBF. CBF 
normally ranges from 3 to 30 Hz. However, 1.5 seconds of data collection time is not sufficient 
time to detect the MWF which is normally in the range 0.2 to. 2 Hz. To analyze CBF and MWF 
for every 1.5 seconds of data collection, two different mathematical analysis tools are needed. 
CBF is analyzed by decomposing the photon count time-series sequence using a wavelet 



BIO1819P0090US 

transfonnation [Strang, G. and T. Nguyen. Wavelets and filter banks. WeUesley-Cambridge 
press. 1996]. Conceptually, wavelet analysis is analogous to Fourier analysis which decomposes 
the signal into sine waves of various frequencies. By replacing the sinusoidal integrand in the 
Fourier transform with a designated waveform of finite duration, better decomposition of the 
fipequencies can be achieved. This waveform is known as a wavelet and the analysis is known as 
a wavelet transfonn. Computationally, this is achieved by digital filter banks composed of levels 
of high-pass and low-pass filters. In digital filters, the filter banks are composed of a series of 
delay operations modulated with coefficients. At each level of iteration, the application of the 
high-pass and low-pass filtering processes to the signal decomposes the signals into the various 
level's of "details" and "^proximations". This process consists of applying a chosen wavelet 
and performing the transfonn at the desired level of iterations (how many levels of high/low pass 
filtering). This results in a set of "detail" and "approximation" coefficients corresponding to the 
hi^ and low pass filter operations. Following the judicious choice of the threshold coefficients 
for filtering, the "detail" signal is dien reconstructed by an inverse transform. CBF is embedded 

A 

in the "detail" signal. 

MWF is analyzed using a Hilbert transformation. In modulation and dranodulation 
communication theory, application of a Hilbert transfonnation to a modulated signal separates 
the modulation signal fixjm the carrio: signal [Schwartz, M. Ihfonhation, transmission and 
modulation and noise: a unified approach to communication systems. 3"* ed. New York: 
McGraw-Hill, 1980. (chapter 4)]. With an assumption that MWF is the modulation signal and 
CBF is the carrier signal embedded in the photon count sequences, the complex envelop will be 
just the low-pass component of ttie analytic signal of the photon count sequences. 
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The approach described herein is superior to existing technology in a number of ways. 
The alignment of the excitation and detection paths using the same microscope objective ensures 
that the excitation focal spot i, the detection focal spot. The use of a long focal length 
microscope objective with long focal depth in conjunction with the conjugate (confocal) spatial 
filter enables the backscattered photons fiom a single focal spot from biological samples such as . 
native epithelia, air-liquid interfece cultures, organ cultures and single cell cultures to be 
detected. The parallel execution of a wavelet transform and a HUbert transform enables the 
respective CBF and MWF to be analyzed simultaneously in real-time from the same group of 
cilia. The present invention envies the interactions and coupling mechanisms between the 
MWF and CBF to be investigated. The system is suitable for a variety of cell motUity 
apphcations, including the beating of cardiac myocytes. 

BRjEEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a drawing of a confocal system for simultaneous measurements of ciUary beat 
frequency and metachronal wave frequency, according to the present invention. 

Figure 2 is a software flow diagram for a photon,' count data acquisition module, 
according to the present invention. . 

Figure 3 is an example of a simulated metachronal ciUaiy activity wave signal composed 
of low and a high frequency signals usmg the modulated LED hght source (LEDPRQtm, 
BioTechPlex). The upper panel is an example of the simulated signal of the ciliary activity. The 
lower panels are the real-time frequencies detennined by the system for two mput frequencies. 
0,4 Hz (3B) and 12 Hz (3C), respectively, using the system. 
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ovine 



Figure 4 is an example of the CBF and MWF detetmined for a sample of native 
tracheal epithelia. The upper panel is an example of the scattered signal of the ciliary activity. 
The lower panels are the real-time determination of the ciUary beat ftequency and metachronal 
wave frequency. 

Figure 5 i§ an example of the CBF and MWP determined for a sample of organ culture of 
ovine tracheal epithelia. The upper panel is an example of the scattered signal of the ciliary 
activity. The lower panels are the real-time determination of the ciliary beat ftequaicy and 
metachronal wave frequency. 

Figure 6 is an example of the CBF and MWP determined for a sample of monolayer 
culture of ovine tracheal epitheUa, The upper panel is an example of the scattered signal of the 
ciUary activity. The lower panels are the real-time determination of the ciUary beat frequency 
and metachronal wave frequency. 

Figure 7 is an example of the CBF and MWP determined for a sample of air-Uquid 
interface culture of ovine epithelia. The upper panel is an example of the scattered signal of the 
ciliary activity. The lower panels are the real-time determination of the ciliary beat frequency 
and metachronal wave frequency. 

The accompanying drawmgs, which are incorporated in and form a part of this 
specification, illustrate embodiments of the invention and, together with descriptions, serve to 
explain the principles of the invention. They are not intended to limit the scope of the invention 
to the embodiments described. It will be ^preciated that various changes and modifications can 
be made without departing from the spirit and scope of the invention as defined in the appended 
claims. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Herein, preferred embodiments of the invention will be described. While the invention is 
described in connection with such embodiments, it will be appreciated that various changes and 
modifications- can be made without departing fi»m the spirit and' scope of tiie invention as 
defined in fbe app&aded claims. 

An embodiment of the optical system for simultaneous measurements, of CBF and MWF 
is schematically shown in Fig. 1. A 0.8 mm beam from a 17 mW He-Ne Laser (Uniphase) [101] 
with hi^ linear stability (99.9% linear) was optically expanded [102] to increase the focal depth. 
An inverted microscope [103] was modified to enable visualization of the beating cilia and 
simultaneous CBF and MWF me^ements. In brief, a 10:90 beam splitter [104] was placed 
along the observation eye-piece optical path [105], A 50:50 beam splitter was placed underneath 
the beam splitter [106]. The expanded laser beam was reflected iq>wards 90 degrees and focused 
by the 40X, long focal length microscope objective (N.A 0.6) [107]. The transmitted laser beam 
was incident on a piece of Teflon [108]. This geometrical configuration enables the 

backscattered photons from the beating cilia which underwent Doppler-shifts together with .the 

I 

photons backscattered from the epithelium which were not Doppler shifted to be collected 
together by th^ microscope objective and further mixed with the non-Doppler-shifted photon 
backscattered from the Teflon. This demodulation process utilizing an additional random 
dqiolarizer to enhance the CBF signal-to-noise ratio is well established in CBF measurements as 
well as in other optical techniques^.. The backscattered light was then reflected by a front-surface 
mirror to a 632 ran laser-light interfer^ce filter (FWHM==10 nm) [109] and spatially filtered by* 
a conjugate pin-hole [110]. The photons were collected with a piano-converging lens onto a 
1.2MHz bandwidth photon countiug photomultiplier tube [ill]. The size of the conjugate 
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pinhole is the size of the primary Airy Disc, S. [Reynolds, G.O, J.B. Develis, G B. Parrent, Jr., 
Brain J. Thompson. Physical Optics Notebook: Tutorials in Fourier Optics. SPIE Optical 
Engineering Press, 1989. Chapter 24]. This is defined as: 
. S = (1.22 X* M *m) / NA, 

where M is the magnification of the objective, m is the magnification of the remaining optical 
components of the system along the optical detection paths; X is the incident laser wavelength; 
N.A. is the numerical aperture of fhe microscope objective. Based on the optical components of 
the converging lens and the beam-splitters, S is estimated to be 25 pm. 

The Doppler shifted photons and the non-Doppler shifted photons collected by the PMT 
yield temporal interf^iice with intensity fluctuations caused by ciliary activity. A PCI bus- 
based, multi-channel A/D; counter timer computer board (PCI 6023E, National Instruments) was 
configured to process .the photon counts sequence for the system [1 12]. Source code was written 
in Matlab, utilizing both the MatLab and C codes. In brief, the photon count events fi'om the 
photodetector associated with the microscope application are sampled as SV transistor-transistor 
logic (TTL) voltage signals by a National Instruments (NT) 6023E counter/timer board residing 
in the PCI slot of the host PC [1 13], 

Referring to Fig 2, the photon count data acquisition module of the software performs the 
following operations: 

a. ) Periodically interrogates the data acquisition module to obtain photon count sequences in 
1024 element arrays. 

b. ) CBF: The raw photon count data in each counting cycle is converted into a sequence of ' 
Os and Is (Clipping) depending on the individual coxmt value being < or > the mean value of the 
raw pjioton count sequence, respectively. 
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Using the MATLAB library function (wavedec), a level 3, 3rd order Daubechies wavelet 
(dbS) decomposition is then performed on the clipped photon count sequence, defined as: 
[C, L] = wavedecpCN.'vmime') 
where: 

X is the clipped photon count sequence 
N = Level of decomposition = 3 
Iwname' = Wavelet type = 'dh3' 
.C = wavelet decomposition output, vector 
L - book keeping ou^ut vector 

MATLAB is avaUable from Mathworks of Natick, M.A. 

Using the MATLAB Kbrary function (wavecoef). a level 1 approximation wavelet 
reconstruction (X') of the cUpped photon count sequence (X) is obtained using the same 
Daubechies (db3) wavelet and the output vectors C and L of the wavelet decomposition process, 
X' is defined as 

X- = wavedec(V,C,L/db3», N') 
where: 

V = Approximation coefficients obtained from the wavelet decomposition procedure 

'db3* = Wavelet used for reconstruction 

N* = Level of reconstruction 

X' is the reconstructed photon count sequence 

Using the MATLAB function Hannitig, a 1024 point symmetric Hanning filter window is 
appUed to the wavelet rw>nstructed photon count sequence, i.e. 
XH' = Conv(Hanning(1024), X*) 
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where: 

Conv = MATLAB Convolution operation 

XH' = Post-Hanning filter, wavelet reconstructed, photon count sequence 

A 1024 point fast Fourier ^transfoim (FFT) of XH' is calculated and the power 
spectrogram of this FFT is obtained. Based on the highest spectral peak component in the power 
spectrogram, the dominant frequency is obtained, classified as ciliary beat firequency (CBF), 
time-stamped and written to a file stored on the computer's hard drive. 

c.) MWF: The metachronal wave frequency calculation uses the same raw photon count 
sequence as the CBF calculation. Firist, the mean photon count value is subtracted from each data 
element of the 1024 element photon count sequence. 

Using the MATLAB fimction Hilbert, a Hilbert transform (XHIL) of the mean subtracted 
photon coxmt sequence (X") is obtained. This is defined as 
XHIL==Hilbert(X") 

A 1024 point FFT is calculated using the absolute value of the ou^ut from Ihe Hilbert 
transform (XHIL). Based on the highest spectral peak component in this power spectrogram, the 
dominant frequency is obtained, classified as metachronal wave frequency (MWF), time- 
stamped and written to a file stored on the computer's hard drive. 

The software also incorporates a graphic user interface (GUI) through which the user can 
control the application and monitor the progress of the data acquisition session including the time 
course of fiie CBF and MWF measur^CTts made during an experiment. The raw photon counts 
are written to one output file and the calculated CBF and MWF values are written to a second 
output file. 
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Example 1; Measurement of simulated metachronal ciliary activity wave signal composed of 
low and a high frequency signals using flie modulated LED light source (LEDPRO*^^* 
BioTechPlex). 

To test the system, a waveform modulated light emitting diode (LED) ligjit source system 
(LEDPRO™, BioTechPlex) was used- A programmable synthesized function gmerator 
(Stanford Res System DS345) was used as an external source to drive a red LED. To simulate 
the metachronal ciliary activity wave signal, a siniisoidal wave fimction selectable from 1 to 40 
Hz, die carrier frequency, is amplitude-modulated witti 80% depth at frequencies ranging fix)m 1 
to .2 Hz. Figure 3a is an example of the raw data obtained with the system. Figure 3b is the 
carrier frequency determined by the system when the modulation frequency is fixed at 0.4 Hz. 

The ^ of the regression line is 0.99. Figure 3c is the modulation frequency determined by the 

. r 

system with the carrier frequency fixed at 15 Hz. The ^ of the regression line is 0.95. 

Example 2: Measurement of the CBF and MWF on native ovine epithelia 

Experimental Protocol 1: Ovine tracheae were obtained from a local slaughterhouse. 
The tracheae were stored in cold Hanks' solution; Each trachea was opened anteriorly. The 
cartilage and the connective tissues were removed by rapid dissection. The resulting posterior 
epithelial tissue of each trachea, comprised of pseudo-stratified columnar ^ithelial cells and 
lamina propria, was stored at 4^C on M199 containmg 0.8% penicillin-streptomycin. The 
epithelial tissues were viable for 5 days. . Figure 4a is an example of the raw photon count data 
obtained with the system. Figure 4b is the CBF determined by the system. Figure 4c is the MWF 
determined by the system. 

1 
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Example s: Measurement of CBF and MWF from organ cxiltured ciliated cells. 

Experimental Protocol 2: The cartilage-free mucosae obtained from experimental 
protocol 1 were cut into small pieces 20 W and washed with Hanks' solution. The mucosae 
were digested at 4°C for 8 hours using 1% protease in M199 containing 1% penicillin- 
streptomycin. Clusters of ciKated cells were harvested by gently shaking the mucosae in M199 
enriched with 10% fetal calf serum and 1% penicillin-streptomycin. The resulting clusters of 
ciliated ceUs were washed with enriched M199 and were added to each collagen coated culture 
chamber. Culture media was changed to Bronchial Epithelial Growth Media (BEGM, 
Clonetics/Biowhitaker). In a humidified, 37^C incubator, the ciliated cells attached to the 
coUagen matrix within one day. The ciliated ceUs displayed columnar morphology and the ciUa 
maintained their beating for longer than 3 weeks. The BEGM solution was exchanged every 3-4 
days. Figure 5a is an scample of the raw photon count data obtained with flie S3«tem (the lower 
panel shows Ihe extended time-scale for MWF). Figure 5b is the CBF determined by the system. 
Figure 5c is the MWF determined by the system. 

Example 4: Measurement of CBF and MWF from single ciliated cells. 

Experimental protocol 3: Resulting clusters of ciliated cells harvested in experimental 
protocol 2 were washed with enriched MI99. They were mechanically dispersed in BEGM 
solution as single ciliated cells. The cell su^ension was added to each collagen coated culture 
chamber. In a humidified, 37^C mcubator, the ciliated cells attached to the collagen matrix 
within one day. The cilia maintained their beating for longer than 3 weeks. The culture solution 
was exchanged every 3-4 days. Figure 6a is an example of flie raw photon count data obtained 
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with the system system (flie lower panel shows the extended time-scale for MWF)... Figure 6b is 
the CBF determined by the system. Figure 6c is the MWF detennined by the system. 

Examples: Measurement of CBF and MWF from a air-liquid interface culture. ' 

Exp^imental protocol 4: Single cells harvested from exp^imental protocol 3 were used. 
About 300 Ml of tracheal /bronchial epithelial cell suspension (1-2x10^ cells) in BEGM with 5% 
fetal calf serum and 1 x lO"' M retinoic acid were dispersed onto a collagen membrane 
comprised of a Cellagen disc (14 mm diameter. ICN Biomedicals). The CeUagen discs are type I 
coUagen penneable membranes mounted on plastic supports. Each of these Cellagen discs 
(culture inserts) was placed ii each of the wells of a 24-weU cell culture plate. 500 ul of BEGM 
was placed on both sides of the discs in each of the wells. The non-adherent cells and cell debris 
of each sample were washed away on the second day. The ceUs were kept in a 37°C, 5% CO2 
incubator with 100% humidity and were allowed to grow in this submerged culture condition. 
Culture medium was renewed every 2-3 days. After five days of submerged culture, cells 
reached confluence and lost their cilia. The cultm« medium of each well was exchanged to 
BECai+lO-'M retinoic acid. An air-liquid interface culture was established by removing the 
medium inside the CeUagen discs and also by lowering the level of medium to the level of the 
Cellagen discs. The cells were kept in a 3rc. 5% CQ2 incubator with 100% humidity. The 
culture medium was renewed every 2-3 days. CiUa re-grew after 7 days in air-liquid interface 
cultures and maintained their beating for the foUowing 2 weeks. Figure 7a is an example of the 
raw photon count data obtained with the system (the lower panel shows the extended time-scale 
for MWF). Figure 7b is the CBF detemiined by the system. Figure 7c is the MWF determined 
by the system. 
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It can therefore be appreciated that a new and novel method and confocal microscope 
system for perfomiing real-time analysis that produces simultaneous tempoml measurements of 
metachronal wave frequency and ciliary beat frequency has been described. While the invention 
has been described in connection with a preferred embodiment tkeceof, it will be appreciated that 
various changes and modifications can be made without dq)arting from the spirit and scope of 
the invention as defined in the appended claims. As a result, the invention is not to be limited to 
the foregoing embodiments, but only to all the appended claims. 
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CLAIMS 

1. A method for simultaneously measuring ciUary beat frequency and metachronal wave 
frequency simultaneously in a sample of native epithelia or cell cultures, such method 
comprising the steps of: 

illuminating said san^jle with an illuminating light beam; 

detecting light backscattered from the sample surface with a microscope system; and 
analyzing the backscattered Kght with a real-time analysis system. 

2. The method of claim 1 in which the illuminating light beam is laser light optically 
expanded to increase its focal depth. ' 

3. The method of claim 1 comprising a further step in which the illuminating Hght beam is 
optically spUt by a beam spKtter such that the transmitted beam is impinged on a piece of 
teflon and the reflected beam is directed along the excitation path of the microscope 
objective. , 



4. 



The method of claim 3 comprising a further step in which the reflected beam produced by 
the beam splitter is spHt by another beam spUtter such that the reflected beam is directed to 
an eye-piece for visuaKzation and the transmitted beam is focused by the microscope 
objective to the sample. ' 
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5. The method of claim 1 comprismg a further step iii which backscattered light from the 
sanq)le and a random scattcrer are mixed. 

6. The method of claim 5 comprising a further step in which the mixed light is collected by a 
photon counting photomidtiplier tube. 

7. The method of claim 1 comprising a further step in which the signals emitted Scorn the 
photomultiplier tube are analyzed to produce simultaneous measurements of metachronal 
wave frequency and ciliary beat frequency. 

8. The method of claim 1 in which the analysis to produce simultaneous measurements of 
metachronal wave frequency and ciKary beat frequency comprises the steps of: 
performing real-time time-scale wavelet transformation of the time series photon count 

data; and 

perfgmiing real-time Hilbert transformation of the time series photon count data. 
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Figure 2 



Phdtodetector 

. — — 

SVTTL Output 
— Slanal 

Counter/Timer 



1024 point photon 
count Sequences 



4 

Photon Count 
CHpping 


i 


Photon Count Mean 
subtraction 


V 




Wavelet 
Decomposition and 
Approximation 




1 

Hilbert transform 






Hanning Window 




Absolute value of 
Hilbert transform 
output 



CBF FFT Power 
Spectrum 



CBF Peak Detection 



Data Storage 



4 1 




Spectrum 




1 r 




MWF Peak Detection 



1 ^ 


Update CBF and 






MWF Data Display 


4 

/ 




20 



Figure 4 
A 



BIO1819P0090US 




21 




22 



I 



BIO1819P0090US 

Figure 6: 
A 




time (sec) 



23 



Figur^ 7: 
A 




BIO1819P0090US 



24 



BIO1819P0090US 

ABSTRACT 

The present invention relates to a method and system that enables continuous real-time 
analysis of both ciliary beat frequency and metachronal wave frequency from a single spot in 
excised native ciliated epithelial tissues as well as in primary and subsequent epithelial cultures. 
Such method and system utilizes the concept of time-scale wavelet analysis and Hilbert 
Transformation for backscattered ligjit derived from a confocal (conjugate) spot on the moving 
cilia. This light contains inherent hi^ and low frequency components corresponding to CBF 
andMWF. 
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